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ABSTRACT 



This report consists of two parts: Part 1, 

Observational study of the time mean flow on transients for 
the northern wintertime for the period 1975/6 to 1982/3, and 
Part 2, Observational study of the divergent flow for the 
tropics and midlatitudes in the northern wintertime for the 
period 1974/5 to 1982/3. The research uses data from the 
tropical global band analysis produced by the Fleet 
Numerical Oceanography Center. 
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Part 1 



Observational study of the time mean flow 
and transients 
for the northern wintertime 

for the period 197-1/75 to 1982/83 
ABSTRACT 

A nine year data set of U. S. Navy operational, twice-daily wind analyses at 200 mb, 700 mb and the sur- 
face from 60°N to 40°S is used to compute statistics on the time mean fields and transient activity during 
the northern winter( December, January, February). The temporal variance and covariance quantities 
were partitioned into low pass (~ 10-90 days) and band pass (~ 2.5-6 days) components. The 
E-vector, as defmed by Hoskins et. al. (1983), and quantities derived from it were the chief diagnostic 
tools. 

The relationships of the transient effects near the exit region of the East Asian and North American 
(200 mb)jet maxima and their associated storm tracks are very similar to the schematic model presented 
by Hoskins (1983), but not for the North African jet maximum. The forcing from'both the synoptic time 
scale and the low frequency eddies is basically equivalent barotropic, with the 200 mb and 700 mb pat- 
terns resembling each other. 

For low frequency eddies, both the subtropical E-vector and the tropical transient kinetic energy 
indicate an equatorward energy propagation and a tropical response at 200 mb over the mid- Pacific where 
the zonal wind is westerly. However, the maxima of these two fields, which are longitudinally aligned, do 
not coincide with the maximum tropical westerlies. An analogous situation is not observed over the 
Atlantic. 

The low frequency calculations reveal barotropic instability regions in the exit regions of the East 
Asian and North American jets and poleward of the NAD maximum. 
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I. INTRODUCTION 



Since the publication of the work of Blackmon et. al. (1977) and Lau (1979 a,b) there has been a 
resurging interest in the time mean aspect of the three dimensional large scale atmospheric flow and vari- 
ability therefrom. Blackmon et. al. (1977) presented data from the National Meteorological Center 
(NMC) operational analyses for 9 winters. The temporal variance and covariance quantities were parti- 
tioned into low pass ( ~ 10-90 days) and band pass ( ~ 2.5-6 days) components. They concluded that 
the jet streams in the time averaged flow develop as a result of thermally direct (time) mean meridional 
circulations over eastern Asia and North America, with thermally indirect circulations over the North 
Pacific and North Atlantic. Lau (1978) using a similar data set found that the transient eddy statistics are 
in good agreement with what would be expected for developing baroclinic waves. The distribution of 
eddy momentum flux is characterized by strong convergence of westerly momentum into the storm track 
axis. 

Both of the above studies and Lau (1979a) have demonstrated that large longitudinal variations 
occur in the time mean and transient eddy quantities. The vertical structure and the transfer properties of 
the bandpass eddies in the storm track areas are typical of developing baroclinic waves whereas those of 
the low pass eddies in the regions of the maximum variance indicate a more equivalent barotropic behav- 
iour. The transient eddy (TE) forcing, which is the largest in the upper troposphere, tends to shift the 
mean subtropical jet (STJ) polewards and aids in counterbalancing the effect of mean vorticity sources 
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and sinks due to surface friction. The main effect of the transient eddy momentum fluxes in the middle 
latitudes is to enhance the barotropic component of the time mean flow and to compensate for the effects 
of surface friction, Holopainen, Rontu and Lau (1982), Holopainen (1983), Holopainen (1978), Holo- 
painen and Oort (198 la, b). 

Recently, Hoskins (1983) and Hoskins et. al. (1983) introduced the concept of an extended Eliassen- 
Palm flux or E-vector as a diagnostic tool to study the transient/mean flow interactions. This vector (or 
more properly pseudovector) provides a measure of the effects of the transient forcing. The traditionally 
used meridional convergence of the zonal momentum flux, <3(u'v')/5y, does not yield a complete descrip* 
tion in the case of zonally varying mean flows. The E-vector contains information (within the limits of 
certain approximations) on the spatial orientation of of the eddies and the direction of their propagation. 
Diagnostic computations of this vector allow a two dimensional description of the effects of transient 
eddies on a horizontal map. Using this vector, Hoskins (1983, 1985), Hoskins et. al. (1983) and James 
(1983) have studied the transient/mean flow interactions in the midlatitudes. 

Studies of the zonally asymmetric time mean circulation must address the maintenance of the sub- 
tropical jets and their interactions with the transient eddies and meridional circulations. However, the 
data coverage in the low latitudes is often sparse and numerical weather prediction (NWP) model forecasts 
are commonly used in the operational data assimilation. Partly because of these reasons previous studies, 
such as Blackmon et. al. (1977), have often excluded tropical regions equatorward of 20°N. Furthermore, 
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all previously published studies using the E-vector computations are limited to data of a single winter sea- 
son. For studies using operational data sets produced by the National Meteorolgical Center (NMC) of 
more than a few years, the analysis system is not homogenous throughout the period. 

In this study we will use the Numerical Variational Analysis (NVA) for the global tropics produced 
operationally by the U. S. Navy to study the time mean and transients of the winter circulations for a 
nine-year period, from 1974/75 to 1982/83. A climatology of the E-vector and quantities computed from 
it will be produced and used as the main diagnostic tools in Part I. For our purpose this data set has 
several advantages, including the absence of any NWP model forecast as a first guess, the retention of the 
analyzed divergence field which has proven useful in synoptic time variation studies over the tropical 
western Pacific, and the stable analysis scheme which has remained unchanged since 1974. The use of the 
tropical data allows us to look at the North African jet in addtion to the attention normally focused on 
the East Asian and North American jets. In addition, it will be possible to compute the equatorward 
energy propagation from the midlatitudes. These were previously studied mainly in theoretical or numer- 
ical models. 

In this paper the data source and computation procedures are described in Section 2. Section 3 
presents the nine-year means of the E-vectors and various quantities computed from it. The results are 
used to discuss the transient mean flow interaction and equatorward energy propagation. Section 4 pres- 
ents the relations of the quantities in Section 3 to the storm track organization and low-level wind shear. 
Section 5 contains concluding remarks and a summary of the results. 
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II. DATA SOURCES AND HANDLING 



A. .FNOC Numerical Variational Analysis 

The wind data used in this work are the operational analyses of the Numerical Variational Analysis 
(NVA) produced by the U. S. Navy's Fleet Numerical Oceanography Center (FNOC). These data are 
produced four times daily by objective procedures on a mercator grid which extends from 40°S to 60°N\ 
The use of the mercator secant projection results in a change in the actual distance between grid points 
from 140 km at 60° N to a maximum value of 280 km at the equator. 

The winds are analyzed at the surface, 700, 400, 250 and 200 mb levels using a successive corrections 
technique based on the work of Cressman (1959). The first guess field is the six hour persistence field. 
These fields are then adjusted to be consistent with a set of numerical variational analysis equations which 
have incorporated the dynamical constraints of the momentum equations with friction included in the 
surface layer, (Lewis and Grayson, 1972). However, the surface and 200 mb wind data serve only as a 
lower and upper boundary condition for the NVA and are not subject to an adjustment. 

An advantage of these data is that no model first guess is used which might bias the analysis in data 
sparse regions. The lack of a model providing a first guess is also a disadvantage since the analysis reverts 
to (6 hour) persistence when faced with a data void. Given the present state of tropical prediction this is 
not an altogether bad situation in the tropics but is considerably less satisfactory in midlatitudes. In any 
case this type of analysis cannot overestimate the transient activity. A small advantage of the Navy oper* 
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ational analyses is that some observations taken from ships are not transmitted on the GTS but are used 
in the analysis at FNOC. Another attractive feature of the data is the fact that the data analysis proce- 
dures have remained virtually unchanged since 1974. Thus the analyses are homogenous with respect to 
the objective analysis technique although the data base has changed as satellite observations have become 
more important and real time aircraft reports have become more numerous. 

III. COMPUTATIONS 

Hoskins et. al. (1983) defines the E-vector as: 



E = -{ (u' 2 - v' 2 ) , (uV) } (1) 



where: 

( )' = deviation from seasonal mean 
( ) = seasonal mean. 

It is important to note that as it stands the quantities used in the above definition are made up of 
transients with periods of 1 to 90 days. In this work a further decompostion of the transients into low 
and high frequency components is made using a low pass filter ( > 10 days ) and a band pass filter (2.5 
to 6 days), respectively. Details of the filters and their response functions are given by Blackmon (1976). 
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Hoskins et. al. (1983) showed that where the E-vectors are divergent (convergent) there is a forcing of 
the mean circulation by the transients with a tendency to increase the westerly (easterly) mean flow. The 
E-vector also contains information relating to the orientation and direction of propagation of the eddies, 
(Hoskins et. al., 1983). A westerly (easterly) E-vector implies meridionally (zonally) elongated transient 
eddies with an eastward (westward) group velocity relative to the mean flow. The E-vector can be con- 
sidered as an effective flux of easterly momentum or the horizontal component of an extended Eliassen- 
Palm flux. It follows that if we define a quantity S m such that: 

V 2 (S m ) = -(VE) (2), 

then in the regions of S m maxima (minima) there is a tendency for westerly (easterly) acceleration of the 
mean flow by the eddies. This quantity rather than the actual divergence is presented, since as pointed 
out by Hoskins et. al. (1983) the E-vector divergence field is noisy and the noise obscures meaningful fea- 
tures. In the calculations the boundary at 60°N does present a problem. This is a region of considerable 
eddy activity, and this boundary will have an adverse effect on the determination of S m The choice of 
zero as a northern and southern boundary when relaxing for S m appeared to give- about as good a results 
as any. 
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Simmons et. al. (1983) showed that the local growth of the kinetic energy associated with deviations 



from the mean state is given by 



d(KE')/dl = CK + CK V 

a y 



where 



CK X = -JL ( O' 3 . l/'-l) I ‘j 

a lcoi(e>) ^ 




CK y - <■'*'( Caifa) J \ + -L_ ivi) 

& L Jo \cos(e) J c°s(e) c)x J 



where u^ and v^ are the components of the time mean wind, a is the radius of the earth, X is the longi- 
tude and 0 is the latitude. 



They further showed that 



d(KE , )/dt - {E*Vu b } 



where u b is the time mean zonal wind. 



For convenience we define a notation such that: 



S E = {E*Vu b } (3) 



In regions where S E is positive (negative) there is a local growth (decay) of the eddies by interaction with 
the time mean flow. 

The transient kinetic energy (TKE) was computed using the formula 



TKE = (u' 2 + v' 2 ) (4) 



This quantity was computed for both the low and band pass transients as well as the unfiltered transients. 
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IV. 200 MB MEAN FLOW / TRANSIENT INTERACTIONS 



A. E-vectors 

Figure 1 presents the mean wintertime 200mb E-vector and zonal wind. In the zonal wind field we 
identify three jet maxima. These will be termed the East Asian jet (EAJ) centered at 135°E, 32°N, the 
North American jet (NAmJ) centered at 75°\V, 35°N, and the North African jet (NAfJ) centered at 
35°E, 27°N. It may be arguable whether the NAfJ should be treated as a separate jet maximum on a par 
with the EAJ and NAmJ or merely as an extended portion of the entrance region of the EAJ. The work 
of Krishnamurti (1961) clearly shows three distinct maxima, for the single winter that he examined. Arkin 
et. al. (1986) have a indication of a third maxima in their five-year average of NMC data. Chang and Lau 
(1980), also based on one season of data, showed that the daily variation of the NAfJ local jet maximum 
is often out of phase with the EAJ. Therefore the NAfJ may be considered as a separate entity from the 
EAJ. On the other hand, examination of individual winter means of the 200 mb wind indicates that in 
some years the NAfJ is rather ill defined compared to the EAJ and NAmJ. 

James and Hoskins (1985) cite evidence that for-the winter of 1980/81 the transient activity about the 
NAfJ is due to internal baroclinic instability. This is in contrast to the EAJ and NAmJ which are regions 
of boundary (Chamey type) baroclinic instability. However, James and Hoskins indicate that the internal 
instability mechanism generates transients that are qualitatively similar with respect to the patterns of 
momentum and heat fluxes to those of the Chamey mode. In the following w'e will attempt to point out 
similarities and differences of the NAfJ system with the other two jet maxima. 
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The E-vectors and transient activity in the Southern Hemisphere are not going to be looked at in any 
detail due to the small amount of data available in that region. This lack of data has an especially detri- 
mental effect on the resolution of the transients when peristence is used as a first guess in the analysis sys- 
tem. This fact combined with the weakness of the summer circulation make any conclusions about the 
Southern Hemisphere transients rather circumspect. 

Figure la is the zonal wind and the low-pass E-vector. As was shown by Hoskins et. al. (1983) for 
the 1979/80 winter, the low-pass transients generally have a westward pointing E-vector .Hoskins et. al. 
(1983) present the 250 mb fields, but we feel that the comparison is valid. This is especially true in the 
exit regions of the East Asian and North American jet maxima and weakly so for the North African sys- 
tem. This implies that these transients are for the most part zonally elongated (u'^ > v' 2 ) and westward 
propagating with respect to the mean flow. The NAmJ and NAl'J are similar in the entrance regions with 
weak northerly E-vectors which implies that v' 2 ~ u' 2 . Figure la may also be compared to the 
E-vectors computed by Simmons et. al. (1983) (their Fig. 16) for the fastest growing mode of their baro- 
tropic model linearized about a climatological mean 300 mb flow. Since the period of this mode is 
approximately fifty days, it may be considered a low pass mode. The two fields are similar in the exit 
regions of the EAJ and NAmJ jets where both show large E-vector magnitudes. However, Figure la 
shows large E-vectors south of the NAmJ which are absent in Simmons' et. al. (1983) mode and Fig- 
ure la does not have the large magnitudes in the Bering Sea region shown by Simmons et. al. (1983). 
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Opposite to the low-pass E-vectors, the band-pass E-vectors ( Figure lb) tend to point eastward. 

2 2 

This implies that the band pass eddies are for the most part meridionally elongated(v' > iT ) and east- 
ward propagating with respect to the mean flow. There are maximum values in the exit regions of the 
EAJ and NAmJ. 

Figure lc is a plot of the unfiltered E-vectors. The direction is predominantly equatorward across the 
jet latitudes. There is a tendency for the low-pass signal to dominate in the EAJ region and for the band- 
pass to dominate in the NAmJ region. The NAmJ and NAfJ are alike in that they have strong northerly 
E-vectors in the entrance regions. In this instance, the NAfJ looks as if it is the entrance region for a 
combined NAfJ - EAJ system. If one takes such a view, the combined system resembles the NAmJ 
structure with respect to the orientation of the E-vectors in the entrance and exit regions. 

The Southern Hemisphere has the reduced transient activity during the summer season. In addition, 
the poor data coverage together with the use of persistence as a first guess undoubtedly underestimates the 
eddy activity there. Nonetheless, in Figure la one can still see the general westward pointing low pass E 
vectors, and in both Figure la and c, the tendency of an equatorward direction. The band pass vectors, 
Figure lb, are too small to make any definite conclusions. 

The E-vectors at 700 mb (not shown) are similar in orientation and relative magnitude to those at 
200 mb. This is most pronounced over the northern portions of the northern hemisphere ocean basins in 
the exit regions of the EAJ and NAmJ. The pattern at 700 mb in the NAfJ region is rather nondescript. 
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The similarity of the 700 mb and 200 mb fields suggests that the 200 mb fields give a good representation 
of the vertically integrated E-vector structure for the northern midlatitudes. 

B. Transient kinetic energy 

Figure 2 presents plots of the low-pass, band-pass, and unfiltered transient kinetic energy (TKE),as 
defined by (4), at 200mb. The gross pattern of the unfiltered field (Figure 2c) resmbles those computed 
by Murakami and L’nninayar (1977, Fig. 3 ) for the winter of 1970-71 and Arkin and Webster (1985, Fig. 
3a) for the average of 1968-79 winters, both based on NMC data. With the exception of the NAfJ which 
has a TKE maximum along its poleward flank, the exit regions of the jet maxima are clearly regions of 
enhanced eddy activity. This could be anticipated from the magnitudes of the E- vectors in Figure 1 For 
all plots the TKE maximum tend to be at about the same latitude of the jet maximum with a tendency 
for a poleward extension. This is consistent with the single season results of Hoskins et. al. (1983), Hos- 
kins (1985) and James (1983). James and Hoskins (1985) noted a maximum of high frequency TKE on 
the poleward flank of the NAfJ between 400 and 200 mb for the winter of 1980/81. They attributed this 
maximum to internal baroclinic instability. 

The model results of Simmons et. al. (1983) consistently point to a low pas? TKE maximum in the 
eastern Pacific, and a secondary maximum in the Atlantic. Both maxima are in the regions of the exit of 
the jets but farther poleward than the maxima of Figure 2a. 
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The low pass and unfiltered TKE (Figure 2a and Figure 2c) display zones of maximum values across 
the Tropics centered at 160°W and around 30°E. This tropical penetration was also noted by Murakami 
and Unninayar (1977). These are both regions of westerlies at 200mb (see Figure 1). Arkin and Webster 
(1985) documented the fact that the equatorial westerlies are a favored region of enhanced TKE. This is 
consistent with the theory of Webster and Holton (1982) which predicted that waves from the mid- 
latitudes could penetrate into the tropics through ducts formed by tropical westerlies. Webster and Hol- 
ton (1982) showed that the observed variance due to low pass transients exhibit maxima in the region of 
time mean westerlies. 

Figure 3 contains plots of low pass TKE, E-vectors and zonal wind for the band 20° N to 20°S over 
the central Pacific. Since the magnitudes of the E-vectors and TKE are generally much smaller in the 
tropics a separate plot of this region is needed to see the details in this region. It also may be used to 
represent the unfiltered pattern because the band pass E-vectors are generally much smaller than the low 
pass in this region in both hemispheres. What is apparent in Figure 3 is the strong equatorward E-vector 
in the central Pacific. While this flux of energy from higher latitudes is coincident with the regions where 
the wind is westerly, the maximum is not in the region of the strongest westerlies. This result is different 
from that of Arkin and Webster (1985), who indicated that the largest TKE values are associated with the 
strongest westerly winds. In Figure 3 the maximum flux from the northern hemisphere is located near 
153W - 160W, which is just to the east of the exit region of the EAJ. The maximum westerlies along 
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15N is, on the other hand, near 130W. Consistent with the maximum E-vectors, the maximum TKE at 



the equator is near 160W. This is a region of westerlies, but again not the strongest westerlies whose 
equatorial location is near 125W. 

The results shown in Figure 3 compare favorably with the modeling work by Hendon and Hartmann 
(1985) and Hendon (1986). Hendon and Hartmann (1985) showed that the maximum variance in regions 
of westerlies in the tropics is due to enhanced equatorward propagation of low pass transients. Hendon's 
(1986) model has only a fairly uniform easterly flow between ION and 10S, and the only forcing is topo- 
graphic. However, this forcing resembles the Tibetan Plateau, therefore the midlatitude and subtropical 
results produced by his model may be considered relevant when compared with the EAJ circulation. His 
model displays strong northerly E-vectors across 20N, just to the east of the exit region of the jet maxi- 
mum at 30N\ As in Figure 3 his E-vectors then turn and become easterly west of 155W and westerly to 
the east of this longitude. 

One possibility that comes to mind in considering Figure 3 is that the tropical transient activity in the 
regions of the tropical westerlies is being driven by transients propagating from the midlatitudes. This is 
probably the case in Hendon's model since it had no convective parameterization. In such a case one 
could sketch out the following relations. The increased(decreased) midlatitude transient activity produces 
stronger(weaker) fluxes (E < 0) equatorward. In regions where the the zonal tropical wind is westerly 
then more(less) energy propagates into the Tropics producing larger(smaller) values of TKE. Figure 4a 
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and Figure 4b are time longitude plots of the seasonally averaged 200 mb unfiltered TKE averaged from 
10°N to 10°S, and the 200 mb Ey( = uV) at 15°N, respectively. There appears to be evidence for the 
association of enhanced tropical TKE and increased northerly E-vector at 160°W. Of course the sample 
size of nine-years renders to speculation how well the atmosphere actually follows the above relation. The 
sources of tropical variability are not driven solely by transients propagating from the northern hemi- 
sphere. In order to evaluate more definitely the amount of tropical variance being forced by midlatitude 
flow a longer time series would be needed so that signal and noise can be separated. 

It is not clear that the same type of midlatitude - tropical interaction is taking place in the exit region 
of the NAmJ near 30°W. Figure 5 shows the low-pass TKE, E-vectors and zonal wind for the tropical 
Atlantic. There is a strong equatorward E-component along the eastern coast of Brazil and a large TKE 
maximum accompanying it. There are northerly E-vectors across the tropical north Atlantic and a TKE 
extrusion along the North Africa coast, coincident with a region of strong westerlies. In any case the pat- 
tern is not nearly as suggestive of an equatorward propagation of transients from the northern midlatitudes 
as was seen in the Pacific, although there is again a strong indication of a flux from the Southern Hemi- 
sphere. 

The difference from the North Pacific pattern is probably a result of the structure of the NAmJ. The 
NAmJ has a northeast-southwest tilt while the EAJ axis is east-west onented. Thus the exit region of the 
NAmJ and the associated TKE activity are farther north. 
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C. Eddy momentum flux 



Simmons et. al. (1983) relates the convergent/divergent patterns of the low pass E-vectors to the sources 

of low frequency barotropic instability. They interpreted strong divergence as an indication of a source 

region of low frequency eddies observed in their model simulations. Figure 6 is a plot of S m { = 

.2 

V (*V«E) } at 200mb. Maxima (minima) of this quantity can be considered as centers of westward 
(eastward) momentum flux convergence due to the actions of the transient eddies. Maxima are regions 
where the transients tend to make the time mean flow more westerly. For the EAJ and the NAmJ the 
effect of the low pass eddies (Figure 6a) is a tendency to decrease the strength of the jet maxima and to 
move the maxima eastward of their observed position. The transient forcing is very weak about the 
NAfJ, however the entrance regions of the NAfJ and the NAmJ look similar. 

The centers of the band pass S m , Figure 6b, are somewhat farther west from the low pass centers for 
the EAJ and NAmJ. They have a tendency to move the jet cores to the north/northeast and to decelerate 
the flow in the jet entrance regions. The band pass S m fields are more zonally elongated and are maxi- 
mum along the poleward side of the jets, therefore they have a tendency to decrease the cyclonic shear of 
the jetstream. The NAfJ pattern is again very weak, but appears to be qualitatively similar to the other 
two jet maxima. 

The unfiltered S m field, Figure 6c, has positive values which extend across both northern hemisphere 
ocean basins. There is a tendency for the EAJ and the NAmJ to have a negative entrance region and a 
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positive exit region. This result is consistent with that of Blackmon et. al. (1977). The NAfJ, on the 
other hand, is more or less surrounded by negative values. 

Hoskins (1983), using scaling analysis, argued that the barotropic forcing should be dominant in the 
jet exit regions. This is borne out in Figure 7 which shows the unfiltered 700 mb S m field. Here the 
maxima coincide with those at 200 mb in the exit regions of the EAJ and NAmJ. 

Many studies of the effects of eddies on the mean flow used the zonally averaged poleward momen- 
tum flux ( u'v' ) only. As shown in Section 3, this is the meridional component of the E-vcctor which 
goes into the computation of S m . We can perform a calculation analogous to S m but just using the ( u 
v' ) component of the E-vector. Denoting the result by the notation S m y, we have 
v 2 ( S my ) = d(-u'v');dy. (5) 

This would produce fields in which the relative maxima indicate the momentum flux convergence due to 
the meridional variations of the (u'v') correlation. The zonal average of the S m and S m y computations 
would be equal. Figure 8 shows the low pass S m y at 200 mb. Compared to the S m shown in Figure 6a 
the S m y maxima in the Pacific are shifted towards North America and generally have somewhat less 
zonal structure than the corresponding S m fields. As reported by Hoskins et. al. (1983), the zonal scale of 
the u'v' correlation is larger than the meridional scale. The zonal variation of u'^ - v'^ becomes largest in 
regions of strong westerlies, leading to large differences between S m y and S m in these regions. Although 
the band pass S m y field (not shown) more closely resembles that for S m , evidently because of less varia- 
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tion of u' 2 - v' 2 in the zonal direction, the unfiltered S m y (not shown) with a rather bland zonal structure 
is quite different from the unfiltered S m ( Figure 6c). These differences point out the importance of zonal 
variations, which cannot be represented by u'v', in considering the local momentum balance. The u'~ - 
v tenn is especially important in determining the local eddy vorticity flux convergence. It has been 
shown in Figure lb that the synoptic-scale transient eddies tend to be elongated in the meridional direc- 
tion (westerly E-vectors). If one is interested in parameterizing the effects of these eddies, a representation 

2 2 

of the time average of the u - v term would be more appropriate than the much discussed uV. 

Another way to examine the transient forcing of the mean flow is to look at the vorticity balance. 
The time averaged vorticity equations can be written approximately as (Hoskins et. al., 1983) 
dQdl + V«VJ; +fV«V + V*{V%') = k«Vx F, 
where f is the Coriolis parameter, and F is the frictional force. 

Taking the vertical integral (denoted by square brackets) through the depth of the troposphere yields 

[d$dt] + ')] ~ [k*Vx fTJ 

Additionally, Hoskins et. al. (1983) show, that to a good degree of approximation: 

V*(V't; ') ^ d(V*E) dy. (6) 

Figure 9 are the low pass, band pass and unfiltered fields at 200 mb of -V'~5(V*E)/<5y. The correspon- 
dence of the low pass and band pass fields with those of Fig. 5 of Hoskins et. al. (1983) is quite good, 
considering that they were looking at only a single winter. 
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The low pass transient forcing is in general greater than that of the band pass. The "quadrapole" 
structure observed in the mid-Paciiic by Hoskins et. al. (1983) is also apparent in Figure 9a and weakly 
discemable in the Atlantic. The vorticity flux convergence is negatively correlated with the mean vorticity, 
indicating that the low pass transients tend to smooth out the time mean (stationary) long waves. The 
band pass forcing Figure 9b, is also similar to that of Hoskins et. al. in that anticyclonic (cyclonic) forcing 
is found equatorward (poleward) of the storm tracks in both northern hemisphere ocean basins. How- 
ever, Hoskins et. al. (1983) found that the band pass vorticity flux forcing is strong only in the Atlantic, 
while Figure 9b, indicates that it is also important in the Pacific. Hoskins et. al. (19S3) note that the year 
they considered was one of anomalously decreased transient activity along the Pacific storm track. 

Hoskins (1983) developed a scaling model for the synoptic transient forcing, in which the vorticity 
flux convergence has a dipole structure that saddles the northern flank of the jet exit regions. The result is 
a tendency for the eddies to move the jet maximum downstream and polewards in the exit region and to 
split the westerly flow at the diffluent end of the jet. The band pass eddy vorticity flux convergence pat- 
terns ( Figure 9b) for the EAJ and NAmJ are consistent with Hoskin's model, but not at all for the NAfJ 
where the net tendency is to move the jet core equatorward. This difference from Hoskin's model is 
apparently due to the smaller zonal extent and, perhaps to a lesser degree, the lower latitude of the NA£J. 

Figure 10 is the mean surface streamlunction computed from the nine year mean surface wind com- 
ponents using the method described in Chang and Lau (1982). Since the vorticity is related to the Lapla- 
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cian of this quantity, it is consistent to think of minima (cyclones) to be sinks of vorticity due to frictional 
processes at the surface and maxima (anticyclones) to be source regions. What is striking is that the eddy 
convergence of vorticity at 200 mb (Figure 9b) is coincident with the surface centers of cyclonic circula- 
tion. This would be expected if the upper tropospheric contribution to the eddy fluxes is the dominant 
term in the vertically averaged, time mean vorticity equation and the the eddy fluxes dominate the verti- 
cally averaged advection. As was shown by Holopainen and Oort (1981b), the balance struck over the 
North Pacific Low is 
V*VV - -g k« V~Xt 

where t is the surface stress. Holopainen and Oort (198 la, b) show that the vertical integral of the eddy 
vorticity flux convergence plays a major role in balancing the surface stress since the advection term tends 
to cancel in the vertical. The synoptic transient eddies act to maintain the surface cyclonic flow to the 
north and anticyclonic flow to the south and east of the storm track. Thus our result verifies the result of 
Holopainen and Oort (1981b) for an independent data set, and also gives a boost in the confidence in the 
transient eddy computations at 200 mb over the Northern Hemisphere oceanic regions. 

D. Barotropic instability (S^) 

Simmons et. al. (1983) integrated a global barotropic model, which is linearized about a climatological 
January 300 mb flow and perturbed by a series of localized forcings in the tropics and subtropics. They 
used a quantity similar to 
Sj: ~ ( E*Vu ), 
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which represents the local conversion of kinetic energy from the mean How to transient eddies, as a meas- 
ure of barotropic instability. On the basis of their results and the work of Hoskins et. al. (1983), they 
suggest that much of the low frequency variability of the Northern Hemisphere wintertime general circu- 
lation is associated with disturbances which derive their energy from the time mean basic state through 
barotropic instability. Figure 1 1 is a plot of Sg for the low pass, band pass and uniiltered data. As noted 
by Simmons et. al. (1983), some caution is needed in interpreting positive values of the conversions as 
indications of barotropic instability. This interpretation is strictly valid only for the global integral. The 
local contribution to the conversion terms do not directly yield local growth of eddy energy because of the 
existence of additional flux terms which distribute the energy while contributing nothing to the global 
mean values. Figure 11a, the low pass Sg, is analogous to Fig. 17 of Simmons et. al. (1983). However, 
important differences are found, which may be expected since the simple structure shown by Simmons et, 
al. is that of a single, most unstable mode. In Figure 11a the large values of Sg in Simmons et. al/s Ber- 
ing Sea region and the North Atlantic near 60N are missing. Instead, the large positive values are con- 
centrated in the belt of 30 - 40N, including the exit regions of the EAJ and NAmJ, and a zone across 
North Africa. Figure 1 la also shows that both the EAJ and the NAmJ have weak negative values in the 
entrance regions. The NAfJ, on the other hand, has a negative Sg in the jet itself. These structures sug- 
gest that the low pass eddies are barotropically generated or maintained in the exit regions of the major 
jets, and north of the NAfJ along 30N. 
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The band pass Sg, Figure lib, values are generally smaller than the low pass. This is not surprising 
as the band pass eddies are known to draw energy from baroclinic instability. In the NAfJ region the 
band pass and low pass fields look similar. On the other hand, the two fields look different elsewhere. 
Positive values are observed in the entrance region of the EAJ and in a zonal belt along 40 N from the jet 
maximum to the exit region. There is no evidence of a secondary positive maximum to the south of the 
jet axis as seen in the model integrations of Hendon (1986). Around the NAmJ, the band pass Sg struc- 
ture is completely reversed from that of the low pass, with North America being a source of kinetic energy 
for the band pass transients and the Atlantic a sink. The unfiltered Sg, Figure 1 lc, reflects the different 
nature of the low and band pass values about the EAJ and the NAmJ, resulting in small values over east- 
ern Asia, the eastern Pacific and most of the Atlantic. 

The different Sg patterns around the three jets, especially those of the band pass eddies, suggest that 
the barotropic instability properties of the three jets are quite different from each other. 

V. STORM TRACKS AND UPPER-LEVEL JETS 

A. Storm tracks 

Figure 12 presents plots of the surface transient kinetic energy for the low pass, band pass and unfil- 
tered wind. The maxima are taken to be the principal storm tracks. The band pass TKE maxima in 
Figure 1 2b have the familiar storm track pattern seen in the data of Blackmon et. al. (1977) and Lau 
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( 1979a, b). It should be noted that Blackmon et. al. (1977) and Lau ( 1979a, b) use maxima in the field of 



1000 mb geopotential variance to define the storm tracks. They occur on the western side of the northern 
hemisphere ocean basins with a southwest to northeast tilt. For both the EAJ and the NAmJ, the band 
pass maximum at the surface is poleward and downstream of the 200 mb jet maximum. However, there 
is a difference between eastern Asia and eastern North America. The TKE contours cut westward deeply 
into the North American continent from the maxima in the Atlantic, but they parallel the Asian coast and 
make the Asian continent relatively storm free. 

The low pass TKE (Figure 12a) in the Pacific has a maxima to the northeast of the band pass maxi- 
mum. This could be explained by the slowing of the storms as they mature and deepen and finally stag- 
nate in the eastern side of the Pacific, forming the Aleutian low. The corresponding feature in the Atlan- 
tic, the Icelandic Low, is off our charts. In the western Atlantic the low pass maxima is almost due north 
of the band pass maximum, this is a common position for storms to stagnate after they have deepend 
quite rapidly while moving northward. The Mediterranean region shows low pass transient activity but is 
quite nondescript in the band pass fields of TKE. 

The uniiltered TKE field (Figure 12c) shows the concentration of activity in the northern part of the 
northern hemisphere ocean basins quite vividly. 

The 700 mb TKE fields (not shown) in general resemble those at the surface. More activity is 
observed in the region of the NAfJ than at the surface. This perhaps in an indication that internal baro- 
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clinic instability may be more important in this region, James and Hoskins (1985). At both levels, there is 
little evidence of any equatorial penetration of TKE in the western Pacific, as was reported by Murakami 
and Unninayer (1977) in their calculation of the 700 mb TKE for the 1970-71 winter. 

B. Low-level wind shear 

Figure 13 is a plot of the magnitude of the low-level (surface to 700 mb) wind shear. To the extent that 
the thermal wind approximation is valid, this is related to the mean low-level thermal gradient. As 
expected, the largest values are found along and to the east of the Asian and North American continents 
beneath the 200 mb wind maxima. There are also large values across the northern African continent. In 
the Southern Hemisphere ( summer ) the strongest gradients tend to be on the west coasts of the conti- 
nents, evidently related to the very cool waters of! these coasts. The positions of the strong shear in the 
western ocean basins is consistent with the interpretation of the band pass TKE maxima in these regions 
as being the result of the growth of baroclinic waves. 

VI. SUMMARY AND CONCLUSIONS. 

Based on the 1974-1983 FNOC global band analysis, a diagnostic study of the effect of transient 
eddies on the winter mean flow was carried out by E-vector and related computations. Our main results 
may be summarized as follows: 
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1. The forcing from both the synoptic time scale and the low frequency eddies is basically 
equivalent barotropic, with the 200 mb and 700 mb patterns resembling each other. 

2. The exit regions of the NAmJ and EAJ are quite similar in structure with respect to many 
of the features due to transients. The entrance region of the NAfJ resembles that of the 
NAmJ, but not the EAJ. In this regard the NAfJ may be considered as the entrance region 
of the combined EAJ/NAfJ system. 

3. The maximum westerly forcing for all transients is on the western end of the storm track to 
the northeast of the major jet (EAJ and NAmJ) maxima at 200 mb. While all transient 
effects are a maximum in the exit regions, the synoptic time scale eddies have a more meri- 
dional orientation while the low frequency eddies are zonally elongated. The difference in 
shapes and the results of the eddy momentum flux convergence indicate that just consider- 
ing the uV correlation will give an incorrect picture of the local momentum balance, and 
that the consideration of the u^ - v'^ term is of vital importance. 

4. The relationships of the transient effects near the exit region of the major jets and stonn 
tracks are very similar to the schematic model presented b? Hoskins (1983), but not for the 
NAfJ. 

5. The eddy vorticity forcing at 200 mb dominates the vertically integrated vorticitv budget. 
This is consistent with the finding of Holopainen and Oort (1983) which suggests that the 
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surface mean circulations are maintained against frictional dissipation by upper level eddy 
fluxes. 

6. For low frequency eddies, both the subtropical E-vector and the tropical transient kinetic 
energy indicate an equatorward energy propagation and a tropical response at 200 mb over 
the mid- Pacific where the zonal wind is westerly. However, the maxima of these two fields, 
which are longitudinally aligned, do not coincide with the maximum tropical westerlies. 
An analogous situation is not observed over the Atlantic. 

7. The low frequency calculations reveal barotropic instability regions in the exit regions of 
the EAJ and NAmJ, and poleward of the NAfJ maximum. Tliis instability distribution 
agrees with Simmons et. al. (1983)'s modeling result in the EAJ, but differs from their 
results considerably elsewhere. 

8. Comparison of the transient kinetic energy fields at the surface indicate the expected tran- 
sition of a synoptic time scale maximum on the western side of the ocean basins to a low 
frequency maximum on the eastern side and farther north. This is consistent with life cycle 
of baroclinic waves developing on the western side and moving north-eastward while deep- 
ening and slowing down. The low frequency transient fields also tend to have an equiva- 
lent barotropic structure over the eastern ocean basins. This is apparently due to the verti- 
cal development accompanied by deepening, (Hoskins, 1983). What is apparent is that the 
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band pass surface TKE is most active in the left exit region of the jet which is to the north 
of the maximum in the low level baroclinicity. The transients have little effect on the 
acceleration of the jet maximum and appear in some regions to have a decelerating effect. 
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Fig. 1A 
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Fig. 1C 
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Fig. 2A 
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Fig. 2B 
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Fig. 2C 
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Fig. 6A 
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Fig. 6B 
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Fig. 6C 
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Fig. 9A 
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Fig. 9C 
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Fig. 12C 
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Part 2 



Observational study of the divergent flow for the tropics and midlatitudes 
in the northern wintertime for the period 1974/75 to 1982/83 

ABSTRACT 

A nine-year data set of U. S. Navy operational, twice-daily wind analyses at 200 mb, 700 mb and the sur- 
face from 60°N to 40°S is used to compute divergent wind during the northern winter(December, Janu- 
ary, February). 

The divergent wind fields at the surface, 700 mb and 200 mb are in good agreement in the Tropics 
and sub-Tropics with the observed OLR fields. The agreement is good in the sense that upper-level div- 
ergence and lower level convergence is associated with low values of OLR, the assumption being that low 
OLR values are areas of vertically developed cumulus. Comparison of the surface and 700 mb patterns of 
velocity potential, indicates that the deepest region of low-level' convergence is found over the maritime 
continent of Indonesia. 

Estimates of the ageostrophic wind using the time averaged vorticity equation at 200 mb are in 
agreement in mid-latitudes with the work of Lau (1979) and Blackmon et. al. (1977). Comparison of the 
ageostrophic estimates and the divergent wmd indicates the former is not a very good estimate of the lat- 
ter. As Blackburn (1985) points out, studies which infer vertical motion from the ageostrophic fields may 
not be on solid ground. 
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I. INTRODUCTION 



Presented here are time mean fields of the divergent wind and velocity potential at 200 mb, 700 mb 
and the surface. These fields can be compared to the results of Arkin et. al. (1986) and Krishnamurti et. 
al. (1973). One aspect here is to compare the analysed divergent circulations to the estimates of ageos- 
trophic wind derived from the time mean zonal momentum budget. The work of Blackburn (1985) indi- 
cates that these two fields should be quite different, and he emphasises the dangers of inferring the diver- 
gent flow from the ageostrophic estimates. Also, it is of interest to investigate how the centers of action in 
the divergent wind field fit into the pattern jet maxima. 

In this paper the data source and computation procedures are described in Section 2. Section 3 
presents the nine-year means of the velocity potential and the divergent wind. The results are compared 
to previous calculations and interpreted in light of what is known of tropical convection. Section 4 pres- 
ents the relations of the divergent centers to the jet maxima and and compares this to estimates of the 
ageostrophic wind. Section 5 contains concluding remarks and a summary of the results. 

II- DATA AND CALCULATIONS. 

The wind data used in this work are the operational analyses of the Numerical Variational Analysis 
(NVA) produced by the U. S. Navy's Fleet Numerical Oceanography Center (FNOC). The details of the 
FNOC analysis procedures are described in Part I. Of importance is the fact that the FNOC procedures 
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are such that no attempt was made to alter the divergence that was analysed in the wind fields. Thus the 
FNOC data represents the longest running operational system containing information about the divergent 
wind. A significant portion of the day to day changes in this divergent wind was found to be consistent 
with physically meaningful processes in work by Chang and Lau (1980, 1982), Lau et. al. (1983) and 
Chang and Lum (1985). 

A. Computation of velocity potential 

The velocity potential ( X ) was computed following the procedures described by Chang and Lau (1980, 
1982). 

In addition, the windfield was directly decomposed into its rotational and divergent components 
using the method of Endlich (1967). This method does not require any assumption about the boundary 
conditions. The divergent wind vectors shown in the following figures are not computed from the velocity 
potential but are those computed by the Endlich technique. The excellent agreement between the X field 
and the divergent winds over the entire grid, gives additional confidence in the accuracy of the computa- 
tions. • 

III. TROPICAL DIVERGENT FLOW AND CONV ECTION. 

Figure 1 presents the nine-winter mean velocity potential and divergent wind at 200 mb, 700 mb and the 
surface. In the tropics the divergence fields are often compared with the outgoing long wave radiation 
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(OLR) under the assumption that low values of OLR imply the presence of deep, persistent convective 
clouds (Gruber and Krueger, 1984). The nine-winter mean OLR is plotted in Figure 2 

The 200 mb divergence has been calculated by several investigators, including Krishnamurti et. al. 
(1973) and Murakami and Unninayar (1977), both for a single winter only, and Arkin et. al. (1986) for a 
five-winter period from 1978/79 to 1982/83. The latter two studies used the NMC operational analysis. 
All these works show a strong east-west (Walker-type) divergent component which is at least comparable 
to the north-south (Hadley-type) component. However, our result shown in Figure 1 indicates a more 
prominent Hadley-type component at 200 mb. This is manifested by a better defined South Pacific Con- 
vergence Zone (SPCZ) as compared with previous studies; the SPCZ outflow corresponds reasonably well 
to the minimum OLR pattern in Figure 2. There are also in general more smaller scale details especially 
compared to the NMC data based calculations. This difference is perhaps due to the influence of the 
NWP model forecasted divergence in the NMC analysis. An example is the outflow zone in the equato- 
rial eastern Pacific centered near 12 1W. This area is totally absent in Murakami and Unninayar (1976) 
and Arkin et. al. (1986), but is present in Krishnamurti et. al. (1973) who based their calculations on a 
subjective analysis of the January -March 1969 data. This divergence feature may^be compared to a nar- 
row band of minimum OLR immediately west of 120W, which is the only stippled ( < 225 Wm“) area in 
the entire tropical eastern Pacific in Figure 2. Other significant differences from the NMC based studies 
include the locations of the divergence center over South America and Africa. They are centered near the 
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equator in the NMC based analyses, while in Figure la they are near 10S, in good agreement with the 
OLR minima in Figure 2 The one feature that shows up in Figure 1 but cannot be identified with an 
obvious corresponding feature in the OLR chart is the smaller scale, secondary divergence center in the 
Arabian Sea/Indian Ocean. 

The more detailed structure in the 200 mb divergence is, in general, consistent with the 700 mb ( 
Figure lb) and the surface ( Figure lc) convergence areas where large-scale deep convection is indicated in 
the OLR. This includes the SPCZ, and the equatorial South American and African convergence zones. 
Among the above mentioned divergence calculations, only Murakami and Unninayar (1976) showed the 
700 mb fields. The main difference shown in our 700 mb results are similar to those found at 200 mb, 
including the more prominent Hadley-type flow, a much better defined SPCZ and greater detail. The 
correspondence of our data with the OLR, especially the 200 mb divergence and the surface convergence, 
is at least an indirect confirmation of the veracity of the divergent wind described in this data. 

There are significant geographical differences in the vertical structure of the low-level divergence. 
Generally, there is a decrease of the divergent w’ind magnitude in going from the surface ( Figure lc) to 
700 mb ( Figure lb). The largest values of £ are found in the equatorial western Pacific. It may be 
worth noting that, in a study of rawindsonde data in summer 1973, Thompson et. al. (1979) found that 
the divergence diminished rapidly with height in the GATE B-scale ship hexagon but does not change 
much in the Kwajalein-Enewetok-Ponape triangle in the western Pacific. It would appear that the two 
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major convective centers over equatorial Africa and South America have their low-level convergence 
largely restricted to below 700 mb. 

An interesting feature of the SPCZ as depicted in Figure 1 is that there is a slightly closer relation 
between the 200 mb divergence and the 700 mb convergence than the surface convergence. Both the 200 
mb and the 700 mb fields clearly show the west-northwest to east-southeast orientation which is not 
obvious in the surface convergence zone. This orientation has been well known from satellite pictures and 
shows up in the OLR. Thus, it appears that certain intrinsic features of the SPCZ are not forced by sur- 
face based processes. 

IV. SUBTROPICAL JETS. 

A. Divergence patterns about the jet maxima. 

Figure 3 is a plot of the 200 mb divergent wind vectors and contours of the 200 mb zonal wind. 
With respect to the upper-level wind features, the EAJ and NAmJ have convergence on the cyclonic side 
of the entrance region and divergence on the cyclonic side of the exit region. The 200 mb and surface 
patterns tend to be of opposite sign indicating deep vertical motions. The vertical motions inferred from 
the 200 mb and surface fields are broadly comparable to the patterns of omega displayed by Lau et. al. 
(1981), This vertical motion field was derived from the operational NMC models which only provided 
data north of 20°N. The NAfJ does not show any of the signs of the divergence patterns associated with 
the EAJ and NAmJ. 
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The divergence maxima in the NAmJ exit is over the TKE surface band pass maxima, Boyle and 
Chang (1987), which is not the case for the EAJ. The TKE surface band pass maximum associated with 
the NAmJ is of greater magnitude than that of the EAJ. The surface convergence zones across the 
northern Pacific and Atlantic are coincident with the surface transient kinetic energy maxima (Fig. 3, 
Boyle and Chang, 1987). This indicates that the convergence associated with the strong cyclones in this 
region dominates in the long term mean. 

In the northern midlatitudes the data of Figure la is in good agreement with the NMC data of Arkin 
et. al. (1986). In both sets of data there are convergence (divergence) centers in the left entrance (exit) 
regions of the EAJ and NAmJ. In Figure la both centers associated with the EAJ lie along 35N, while 
the centers of the NAmJ tilt southwest to northeast from 35N to about 52N. This is in keeping with the 
orientation of the jet axis. The velocity potential and divergent winds of Lau(1979, Fig. 13) is similar to 
Figure la for the NAmJ but his data show strong northerly divergent wind in the EAJ exit, while our data 
( and Arkin et. al., 1986) have easterly to southeasterly divergent winds.' Lau (1979) has a major conver- 
gence center at 170E, 25N but is missing the divergence center east of the dateline at 35N seen in Arkin et. 
al. (1986) and in Figure la. Lau's (1979) values were deduced using a vorticity budget technique on the 
NMC wind analysis. Lau could only compute values poleward of 25N due to data limitations. As men- 
tioned in the section above the fields of Arkin et. al. (1986) and Figure la. are quite different in the 
Tropics, yet they are quite similar in midlatitudes especially in the disposition of the convergence/ 
divergence centers about the major jet maxima. 
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B. Divergent versus ageostrophic flow. 



The. acceleration of the zonal wind may be expressed as 



D(u)/Dt = fv a 



or 

du /dt + xxdu/dx + vduidy ~~ f(v - v ). (1) 

For a long term average, duidt is approximately zero. The actual balance as shown by Lau (1979) and 
Blackmon et. al. (1977) is mainly between the zonal advection and the coriolis torque. 

Given that the flow about the jets is for the most part horizontal and zonal at 200 mb and that the 
divergent wind is a good approximation to the ageostrophic wind, one would expect poleward divergent 
tlow in the entrance regions to the jet maxima (where udu/dx < 0 ) and equatorial divergent flow in the 
exit regions. That this is generally the case is apparent from Figure la, and Figure 2. The EAJ does 
appear to have southerly divergent flow in the entrance region but no clearly organized northerly values in 
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the exit region. The NAfJ has evidence of positive values in the entrance and negative in the exit. The 
NAmJ has definitely positive entrance values if the entrance is considered as extending to the southwest of 
the Baja Peninsula. 

Of course, expecting such a simple picture about the time mean jet is rather naive. Blackbum(1985) 
has shown a significant amount of the ageostrophic flow is contained in the rotational component of the 
wind. He has shown that for the winter of 1982, 83 the cross contour flow which forms a couplet about 
the EAJ is for the most part is due to the rotational wind. That the rotational wind is important in the 
local generation of kinetic energy was shown by Krishnamurti (1968). 

The time mean ageostrophic flow can be estimated using the zonal momentum equation as follows: 

f(v ■ v a) " u<5u/<3x + vdu/dy + d(u'u')dx + d(u'v')dy (2) 

© 

where the overbars indicate a long term time average and du/dl is assumed to be negligable. The right 
hand side of (x) was computed and the results are displayed in Figure 4. The pattern is very similar to 
that of Lau (1979) who performed a similar computation on a different data set. The fields about the jet 
maxima are dominated by the udu dx term, as found by Lau. For comparison Figure 5 presents the 
product of the corioiis parameter and the meridional divergent wind (fv^). The v^ was computed by the 
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method of Endlich described in Section 3. There is some similarity between the fv a and fv^ fields in the 
subtropics but the correspondence degenerates quickly poleward. Even in the subtropics the agreement 
across Northern Africa is non-existent. This poor correspondence between the fields could be anticipated 
from the results of Blackburn (1985). This highlights the dangers of inferring the divergent (vertical 
motion) circulations from estimates of the ageostrophic wind. 

V. SUMMARY AND CONCLUSIONS. 

1. The divergent wind fields at the surface, 700 mb and 200 mb are in good agreement in the 
Tropics and sub-Tropics with the observed OLR fields. The agreement is good in the 
sense that upper-level divergence and lower-level convergence is associated with low values 
of OLR, the assumption being that low OLR values are areas of vertically developed 
cumulus. Comparison of the surface and 700 mb patterns of velocity potential, indicates 
that the deepest region of low-level convergence is found over the maritime continent of 
Indonesia. 

2. Estimates of the ageostrophic wind using the time averaged vorticity equation at 200 mb 
are in agreement in mid-latitudes with the work of Lau (1979) and Blackmon et. al. (1977). 
Comparison of the ageostrophic estimates and the divergent wind indicates the former is 
not a very good estimate of the latter. As Blackburn (1985) points out, studies which infer 
vertical motion from the ageostrophic fields may not be on solid ground. 
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